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COARSE-GRAINED TRANSPORT OF A TURBULENT FLOW VIA MOMENTS 
OF THE REYNOLDS-AVERAGED BOLTZMANN EQUATION 


RAFAIL V. ABRAMOV 

Abstract. Here we introduce new coarse-grained variables for a turbulent flow in the 
form of moments of its Re 5 molds-averaged Boltzmann equation. With the exception of the 
collision moments, the transport equations for the new variables are identical to the usual 
moment equations, and thus naturally lend themselves to the variety of already existing 
closure methods. Under the anelastic turbulence approximation, we derive equations for 
the Re 5 molds-averaged turbulent fluctuations around the coarse-grained state. We show 
that the global relative entropy of the coarse-grained state is bounded from above by the 
Re 5 molds average of the fine-grained global relative entropy, and thus obeys the time 
decay bound of Desvillettes and Villani. This is similar to what is observed in the rarefied 
gas d 5 mamics, which makes the Grad moment closure a good candidate for truncating the 
hierarchy of the coarse-grained moment equations. We also show that, under additional 
assumptions on the form of the coarse-grained collision terms, one arrives at the Navier- 
Stokes closure, which can be naturally extended to the Burnett and super-Burnett orders. 
Finally, we suggest crude parameterizations of the coarse-grained collision terms for use 
as starting points in numerical simulation and modeling. 


1. Introduction 


A common principle in turbulence modeling involves the averaging of the well-known 
Navier-Stokes equations to filter out the rapid small-scale turbulent fluctuations from 
their solutions to reduce computational cost. This is usually done via the Reynolds aver¬ 
aging of the Navier-Stokes equations which produces additional nonlinear 

terms in the equations. These additional terms are often modeled via the turbulent eddy 
viscosity assumption, which was suggested back in 1877 by Boussinesq |j^, and which 
generally provides good agreement with experiments in near-wall boundary layers. 

However, for the fully developed three-dimensional turbulence it was observed that 


the concept of the turbulent eddy viscosity fails |19|. Girimaji [[^| argued fhaf the 
reason why the turbulent viscosity approximation could be invalid for the averaged 
Navier-Stokes equations was that the latter required that the corresponding solution 
of fhe Bolfzmann equation |[^|7 13 21| was near ifs Maxwellian equilibrium, while the 


corresponding averaged solution did not have to be near equilibrium. As an alternative 


to the Reynolds averaging of the Navier-Stokes equations, Girimaji |12| proposed the 


direct filtering of the solution of the Boltzmaim equation instead, and then solve the 
filtered equation using the lattice Boltzmaim method (LBM). 
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However, Girimaji's approach had a few drawbacks. First, the Reynolds stress closure 
problem was not eliminated in p^ , and the Reynolds stress was modeled by the stan¬ 
dard Smagorinsky-Lilly closure |33| , the assumptions of which can be traced back to the 
Boussinesq approximation Second, the Bhatnagar-Gross-Krook (BGK) approxima¬ 
tion of the collision terms Q, which was used in |[T^, sets the Prandtl number of the 
flow strictly to 1, whereas it is known that the Prandtl number of a fluid is generally dif¬ 
ferent from 1; for example, for a monatomic ideal gas it equals 2/3, and its value for the 
air is around 0.7-0.8 (there are, however, improved BGK collision parameterizations with 
non-unitary Prandtl numbers lj^|^)- Third, the lattice Boltzmaim method usually in¬ 
volves more computational variables than the standard fluid dynamics methods, which 
may limit its use in some applications. 

In this work we propose a new coarse-graining approach were the Reynolds averag¬ 
ing is used directly on the Boltzmaim equation like in p^ , however, we further convert 
the resulting Reynolds-averaged Boltzmann equation into the hierarchy of the coarse¬ 
grained moment equations. While the new coarse-grained moments are different from 
the usual Reynolds averages of moment variables used in conventional transport meth¬ 
ods, the resulting hierarchy of the transport equations for the coarse-grained moments 
is identical to the usual moment transport hierarchy, with the exception of the nonlin¬ 
ear collision terms. We show that the global relative entropy of the Reynolds-averaged 
Boltzmaim equation is bounded from above by the Reynolds average of the global rel¬ 
ative entropy of the usual Boltzmann equation, and, therefore, obeys the same bound 
on the decay rate as established by Desvillettes and Villani p^ for the usual Boltzmann 
equation. This justifies the same closure methodologies as are used to truncate the mo¬ 
ment transport equations for the rarefied gas dynamics 0, in particular the Grad p4p^ 
closure. 

For the modeling of statistical properties of the turbulent fluctuations, we derive the 
equation for the transport of the turbulent kinetic energy under the approximation of the 
anelastic turbulence, as well as the appropriate relations for the turbulent components of 
the stress and heat flux. There remain three unspecified dissipation terms: the turbulent 
energy dissipation rate, and the two collision moments for the coarse-grained stress 
and heat flux. The reason why the coarse-grained collision terms remain unspecified 
is because they originate from the Reynolds average of the nonlinear collision operator 
of the Boltzmann equation, and thus require additional information about the structure 
of the statistical ensemble (and, therefore, the physics of the flow). We suggest crude 
approximations for these remaining nonlinear terms, based on a dimensional reasoning. 

The manuscript is organized as follows. In Section]^ we discuss the Boltzmann equa¬ 
tion and its conventional closures: the Euler, Navier-Stokes and Grad equations. In 
Section we introduce the Reynolds averaging operator, define the new coarse-grained 
variables, and show that the transport equations for the new coarse-grained variables 
have the same hierarchy as those for the conventional fine-grained variables, with differ¬ 
ent collision terms. In Section]^ we parameterize the turbulent fluctuations of the stress 
and heat flux, and derive the transport equation for the turbulent energy. In Section 
we describe the Grad and Navier-Stokes moment closures for the new coarse-grained 
transport equations, as well as mention how to extend the latter to the Burnett and 
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super-Burnett orders. In Section we discuss basic approaches to the modeling of the 
coarse-grained collision terms and the turbulent energy dissipation rate. In Section]^ we 
summarize the results and discuss future work. 


2. Conventional moment closures of the Boltzmann equation 


In the absence of external forces, the Boltzmaim equation for a 3-dimensional fluid is 
given by 

(2.1) '^+vVJ = C{f). 

Here, t and x are the time and space coordinates, v is the velocity of a fluid particle, 
and f{t,x,v) is the statistical velocity distribution of the fluid particles, at the location x 
and time t. The left-hand side contains the transport terms for /, while the right-hand 
side contains the nonlinear collision term, which generally has the effect of dissipation. 
In many applications, C (/) is assumed to be bilinear in /, as the situations where three 
or more particles collide at once rarely occur. 

Let the angle brackets denote the average over the fluid particles v: 

(2.2a) {b)f{t,x) = J b{v)f{t,x,v) dv, 

(2.2b) {b)c{f)it,x) = J b{v)C{f{t,x,v))dv. 

for an integrable function b(v). Then, one can apply these averages onto the Boltzmarm 
equation in ( |2.1| > and obtain the corresponding moment transport equation: 

g 

(2.3) -{b)f + div^{bv)f={b)c(f). 

The nonlinear collision term C{f) has the requirement of the mass, momentum, and 
energy conservation: 

(2-4) (l)c(/) = 0, (u)c(/) = 0, {\\ vf ) cif ) = 0. 


The last identity also signifies that there is no "internal energy" in the fluid - all energy 
that the fluid carries is confined to the velocity of its particles. This transport-collision 
model applies, for example, to a monatomic ideal gas. Here we adopt this model for 
the sake of simplicity of illustration, as the presence of "hidden" energy-accumulating 
degrees of freedom in the fluid particles requires a different, more complicated treatment 
of the Boltzmarm equation, which will be presented elsewhere. 

In order to obtain the conventional moment closures from the moment transport equa¬ 
tion (2.3|>, we first introduce the following conventional velocity moments of /: 


(2.5a) 

p = (1) f, density. 

(2.5b) 

pu = {v)f, momentum. 

(2.5c) 

p = ^( u —M pressure. 
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(2.5d) 

(2.5e) 


S = {{v — u) ® {v — u))f — pi, stress, 
| 2 / 


q =-{\\v — u\f{v — u))f, heat flux, 

where denotes the outer product of two vectors. Rearranging the averages ( 2.5c| >- 
( |2.5e| > above, one writes the identities 

(2.6a) 


\v\ 


(2.6b) 


(2.6c) 


= 3p + p\\uf, 

{v ® v) f = S + pi + pu ® u, 


1 




'v) f = q + Su+ xpw + -| 0 ||m|pm. 


Now, writing the moment transport equation in (2.31 for the velocity moments in (2.Sal, 
( 2.5b| |, ( 2.6a H 2^^ / expressing those moments in terms of p, u, p, S and q using the 
identities above, and taking into account the mass, momentum, and energy conservation 
laws in (2.4||, one arrives at the following system of equations: 


(2.7a) 

(2.7b) 

(2.7c) 

(2.7d) 

as 


dt 


^ + div(|OM) = 0, 


+ div(| 0 M u + pi + S) =0, 


^ + div(pM) + ^(pdivM + S : (V 0 m) + div^) = 0, 


— + (m • V)S + div(M)S + 2[S(V 0 m)] jq + divQ - -(div^)I + 2p[V = Cs, 


dt 

(2.7e) 

dq 

dt 

where 


^ + div(^ 0 m) + (^ • V)m - - ( ypf + S^ div(p7 + S) + Q : (V 0 m) + divR 


( 2 . 8 ) 




Cs = {V®v)c^f), Cc! = -{{\\vf -2u-v)v)c^fy 


— Cfl, 


the symbol denotes the Frobenius product of two matrices, and [A] denotes the 
traceless symmetrization of a 3 x 3 matrix A: 

[A],, = t(A + Aq-tTr(A)I. 


(2-9) TS 2 . 

Above in ( 2.7| |, Q and R are the unknown higher-order moments, 
(2.10a) Q = {{v — u) ® {v — u) ® {v — u)) f 

being the full 3-rank skewness tensor, and 


(2.10b) 


R= -{\\v — uW'^{v — u) ® {v — u))f 
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being the matrix of the contracted 4th-order moment. Both Q and R obey their own 
transport equations, which, in turn, obviously include moments of yet higher orders, 
and so forth. 

In order to understand how to close the moment equations above, we turn to the well- 
known Boltzmaim's H-theorem for ideal gases. We state Boltzmaim's H-theorem in the 
same form as in Golse 


Proposition 1 (Boltzmaim's H-theorem). The following inequality holds for the collision term: 

(2.11) (ln/)c(/) < 0. 

Moreover, the following three conditions are equivalent: 

(ln/)c(/) = 0, 

C{f) = 0 for all V G R^, 

/ is the local Maxwellian distribution, 

where 9 denotes the temperature 

e=P. 

p 

In order to make use of this theorem, one introduces the local and global entropy 
functionals as follows. The local entropy Si [f] (f, x) is given by the functional 

(2.14) Si[f]{t,x) = - J flnfdv = -{lnf)f. 

The global entropy Sg[f] (f) is further given by 

(2.15) Sg[f]{t) = J Si[f]dx=- J {\nf)fdx. 

Now, we can look at the evolution equation for (In/) y, which is given by 

^ /(ln/ + l) {-V -Vxf + C{f)) dv. 
Observing the mass conservation law in ( |2.4| >, we further derive 

(2.17) ^(ln/)^ + div;,(z;ln/)/= (ln/)c(/). 

Assuming that one can omit the effect of the spatial boundaries when integrating, we 
obtain for the global entropy 

(2.18) = - J (ln/)c(/) dx. 


( 1 ) 

( 2 ) 

( 3 ) 

( 2 . 12 ) 


(2.13) 


Qualitatively, the equations in ( |2.17[ ) and ( |2.18| > do, roughly, the following: 
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(1) The equation for the local entropy in ( 2.17| > tends to increase the local entropy 
S; [/], unless / is already the local Maxwellian state in ( 2.12| >, in which case S/ [/] is 
already at its maximum for given local constraints p, u and 6 (or, equivalently, p). 
The rate of convergence towards the local maximum entropy state here is usually 
very rapid, however, the presence of advection prevents S/[/] from reaching its 
local maximum at (|2.12|). 


(2) The equation for the global entropy in ( |2.18| > tends to increase Sg[f], unless / is 
the global Maxwellian distribution of the form 


(2.19) 


/m(^) = 


Po 


( 27100 ) 3/2 


exp 


\V - Mol 


200 


where pQ, uq and 0o are constants throughout the spatial domain, specified by 
the total mass, momentum and energy constraints in the system. This process is 
unaffected by advection, and the time rate of convergence to the global maximum 
state is (Desvillettes and Villani |10|), that is, rather slow. 

These two processes are, in a certain sense, "mutually exclusive", that is, the closer the 
S[f] {t,x) is to its local maximums at each x, the slower the rate with which Sg[f] (f) ap¬ 
proaches its global maximum. This relation between the two processes can be quantified 
more systematically by introducing the relative entropy (Kullback-Leibler distance, p0| ) 
H[f\g] between two distributions / and g as 


( 2 . 20 ) 


H[/lxl(0 = JfMf/g) 


dv dx. 


It is easy to show that H[f\g] is always non-negative, and is zero if and only ii f = g. 
Now we denote the local and global relative entropies as 


( 2 . 21 ) 


H||/I=H[/I/m1, Hg\f] = H [f Iff,], 


respectively. Then, Hi [f] and Hg [f] measure the information-theoretic distance between 
/ and ( |2.12[ ) or ( |2.19| >, respectively, becoming quantitative indicators of the relation be¬ 
tween the local and global entropy processes, described above. The "mutually exclusive" 
behavior (slow decrease rate of Hg for small values of Hi, and, vice-versa, rapid decrease 


of Hg for large values of Hi) can be observed in Figure 5 of |10|, and also in Figures 8 


and 9 of Filbet, Mouhot and Pareschi |111, where the time series of the local and global 
Kullback-Leibler distance between / and both ( 2.12| > and ( 2.19[> a re displayed for a direct 
numerical simulation with a Boltzmaim equation. In Figure we show an example of 
such time series, adapted from Figure 8 of |11|. In fact, if, in a hypothetical situation, 
Si[f] is always at its maximum at all points x (which amounts to strictly zero H;[/]), 
then Sg[f] must be constant in time even if it is below its global possible maximum (or, 
equivalently, Hg[f] must be constant in time even if nonzero). This assumption is, of 
course, unrealistic in practice, however, it is used in the Euler closure of the moment 
equations as we show below. 

Different assumptions about the time-scale separation between the local and global 
entropy processes lead to three different conventional moment closures of the moment 
transport equations in ( |2.7| >. 
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Figure 1. Time series of the local and global relative entropies of a solution 
of the Boltzmann equation. This plot is adapted from Figure 8, upper- 
right plot, of the article "Solving the Boltzmann equation in N log 2 N" by 
F Filbet, C. Mouhot, and L. Pareschi, published in SIAM J. Sci. Comp., 
28(3):1029-1053, 2006 (see reference [111 in the bibliography section). Used 
with permission of F. Filbet, C. Mouhot, L. Pareschi and the Society for 
Industrial and Applied Mathematics. 


2.1. The Euler closure. In the Euler closure, it is assumed that time-scale separation 
between the local and global entropy processes is infinite. Namely, it is assumed that the 
local entropy S;[/] always instantaneously jumps to its maximum for given p, u and 6, 
and f{t,x,v) is permanently the Maxwellian state in ( 2.12| >. The global entropy Sg [f] is, 
therefore, fixed at its initial value, irrespective of its maximum state for the present total 
mass, momentum and energy constraints. This assumption allows to set the stress S and 
heat flux q identically to zeros everywhere, since their corresponding velocity moments 
are indeed zeros for the Maxwellian distribution in ( |2.12| >. The resulting famous Euler 
equations are given by 


(2.22a) 


^ + div(| 0 M) = 0, 


(2.22b) div(pu C) u) + Vp = 0, 

(2.22c) ^ -|- div(pM) -|- divM = 0. 

The practical application of the Euler equations is somewhat limited due to the lack of 
any kind of diffusion. 
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2.2. The Navier-Stokes closure. The Navier-Stokes closure is more sophisticated than 
the Euler closure in the sense that it recognizes the finiteness of the time scale of the 
local entropy dissipation process, as well as takes into account the correction in the local 
entropy process due to advection. As a result, dissipation of the global entropy is present 
on the long time scale, which is more physically realistic than the constant global entropy 
in the Euler closure. 

As a result of statistical analysis of collisions of large numbers of ideal spheres (which 
ideal gas is modeled upon), one arrives at the following first-order linear approximations 
of the stress and heat flux collision operators in the equations for the stress ( 2.7d| > and 
heat flux dZTeb: 


(2.23) 


Cs ^ -^S, 

}i 


5 p 


where the coefficients p. and k are the viscosity and heat conductivity of the fluid, respec¬ 


tively (for more details, see, for example, |Sl4|). Both y and k are defined largely by 
the physical properties of colliding spheres, and for ideal gases they also weakly depend 
on the temperature 6. The ratio of the heat flux over stress dissipation rates is known as 
the Prandtl number. 


(2.24) 


Pr = 


heat flux dissipation rate 5 y 
stress dissipation rate 2 k ' 


which for monatomic gases equals precisely 2/3 

In what follows, we prefer to keep both the viscosity y and heat conductivity k as 
unrelated parameters, in order to distinguish between the terms related to the stress, 
and those related to the heat flux. However, it must be noted that, strictly put, all 
derivations below are only valid when 15^ = 4 k, since this is what the monatomic gas 
model inherently implies. 

Both y and k are proportional to the masses of colliding spheres, which, in the case of 
molecules, renders both y and k very small. As a result, the time scale of local entropy 
dissipation is ^ y/p (or, equivalently, ^ K/p), and if p is not too small, both the stress 
S and heat flux q are very rapidly damped by their collision operators towards their 
Maxwellian values (that is, zeros). This allows to simplify the transport equations for 
the stress ( |2.7d[ ) and heat flux ( |2.7e| > into 


(2.25a) 


9S _ r_ . P „ 

3^+2p[V®«],,«--S, 


(2.25b) 


dq 

dt 





5 p 
2k 




where the advection terms with S, q and Q were dropped, and the contracted fourth- 
order moment R was replaced with its Maxwellian value 


(2.26) 


R 


5p^j 

2p 
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Now, assuming that ^ and k are small enough so that S and q decay on a much shorter 
time scale than the evolution of p, u and p, one then approximates the stress and heat 
flux by their approximate steady states: 


S = —2p [V (8) u] 


TS' 


q = -KV{p/p). 


(2.27a) 

(2.27b) 

A similar procedure is described by Grad pA| , pp. 371-372, and by Struchtrup and 
Torrilhon | |^ . 

The important consequence here is that the attracting states for S and q are not ze¬ 
ros (as opposed to the Euler closure) due to the effect of advection, which was taken 
into account. Lastly, substitution of the approximate states in ( |2.27[ ) into the transport 
equations ( |2.7a[ )-( [Z^ yields the famous Navier-Stokes equations: 


(2.28a) 


(2.28b) 


(2.28c) 


dt 


^ + div(|OM) = 0, 


div(|0M 0u) + Vp = 2 div (^ [V G) m] 


+ div(pM) + -p dlVM = -p 


[V G) u] 


TS 


I TS 


+ ^div(KV(p/p)). 


Observe that now the Laplace diffusion is present in both the momentum and pressure 
equations, on the characteristic time scale of ~ (or, equivalently, ~ R~^)- This 
Laplace dissipation has the effect of "smoothing out" the velocity and pressure fields, 
which tends to slowly increase the global entropy Sg[f] (or, alternatively, decrease the 
global relative entropy Hg [f ]), making the process more consistent with the findings of 
Desvillettes and Villani |10| . The term with p in the pressure equation is "ill-posed" (in 
the sense that it tends to increase the pressure, without any counter-balancing), and thus 
is often dropped. 

In what follows, we assume that the Navier-Stokes relations ( 2.27| > are valid approxi¬ 
mations for the fine-grained stress S and heat flux q at all times. 

2.3. The Grad closure. Under certain conditions, the pressure p becomes so small that 
the characteristic time scales ^ p/p and ^ p“^ become comparable. This happens 
when the fluid is "rarefied", that is, the fluid particles are so spread out in space that 
the collisions rarely happen 0. In this situation, the local entropy S/[/] and the global 
entropy Sg[f] may exhibit comparable time scales of evolution, and the Navier-Stokes 
(let alone Euler's) closure becomes inapplicable. In particular, one can no longer assume 
that the solution / of the Boltzmaim equation in (2.1) is a local Maxwellian ( |2.12| ), or 
anywhere near it. 


An appropriate closure for this situation was suggested by Grad |14,15|. The Grad 


approximation is based on the Hilbert expansion of the distribution function / in the 
Boltzmarm equation around its local Maxwellian state (2.12), such that it has the pre¬ 
scribed stress and heat flux, in addition to the density, momentum and pressure. The 
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resulting distribution approximates the statistical state of the fluid away from the local 

thermodynamic equilibrium and is given by 

(2.29) 


fciv) = fM{v) 


1 + ^ 


\v — u\ 


5p/p 


-l]q-{v-u) + : ((u - u) ® {v - u)) 


where is the corresponding Maxwellian distribution ( |2.12| >. As a result, the full skew¬ 
ness tensor Q and the contracted flatness matrix R are approximated by their values 
provided by ( |2.29| >, which can be computed explicitly as: 


(2.30a) 

(2.30b) 

In particular, this means 
(2.31a) 






divQ = ^ [V 0 + ^(div^)I, 


(2.31b) 


Q'.iy®u) = - ((divw)^ + (^ • V)m + (V 0 u)q). 
5 


Substituting these Grad approximations into the transport equations ( |2.7| >, we obtain the 
full set of closed transport equations for the coarse-grained variables p, u, p, S, q: 


(2.32a) 

(2.32b) 

(2.32c) 


dt 


^ + div(pM) = 0, 
div(pM 0 M pi S) =0, 


^ div(pM) + ^ (P divM S : (V 0 m) div^) = 0, 


as 


(2.32d) (m • V)S div(M)S 2[S(V 0 m)] - [V 0 2p[V = Cs, 

(2.32e) ^ ^(div(^0M) (q ■ V)u) + ^((V ®u)q - {u • V)q)- 


1 f5 


pi + S^ div(p7 + S) + div (^^1 






Observe that the collision terms Cs and Cq retain their general form, while before in the 
Navier-Stokes closure they were approximated by the linear damping terms. Although 
the equations in ( 2.32| > are known to develop shock waves at high Mach numbers | T4}T5|, 
a suitable regularization have already been developed by Struchtrup and Torrilhon |34 
35 by applying a similar type of closure as in ( |2.25| >-( |Z^ , but to the higher-order 
transport equations for Q and R, which produced additional diffusion terms in the 
equations for the stress and heat flux. 
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2.4. Higher-order moment closures. Obviously, one does not have to close the hierarchy 
in ( |2.7| > at the equations for the stress S and heat flux q, and instead choose to include 
the transport equations for Q, R and higher-order moments However, 

for the clarity of presentation, we restrict ourselves here to the relatively simple closures 
presented above. 


3. The moments of the Reynolds-averaged Boltzmann equation 


In a turbulent flow, the quantities p, u, p, etc usually exhibit rapid small-scale os¬ 
cillations in both space and time (that is, they are "noisy"), which often renders the 
direct numerical simulation of such flow computationally intractable, due to the need 
of extremely fine space and time discretization. Instead, one resorts to solving a differ¬ 
ent system of transport equations, which resolves the large scale features of the flow while 
forgoing the small scale oscillations in favor of a more smooth (and, of course, approx¬ 
imate) solution, in order to afford a more coarse computational discretization. Below, 
we refer to such approximate "smoothed out" variables as the coarse-grained variables, 
superscripted by asterisks. 

Historically, such coarse-grained flow equations where obtained by applying the oper¬ 
ation of "averaging" onto the turbulent Navier-Stokes equations, which was purported 
to smooth out small scale rapid oscillations. Reynolds pioneered the use of such av¬ 
eraging operators back in 1895, and they became therefore named after him. Throughout 
the twentieth century, multiple different interpretations of the Reynolds averaging oper¬ 
ator were developed, including time averaging tor steady coarse-grained flows, running 
time-average for unsteady flows, scale-dependent spatial averaging (which later evolved 
into the large eddy simulation method |j^, or LES), and so forth. Eventually, the idea 
of statistical Reynolds averaging was adopted (see |^26| and references therein), due 
to its convenient mathematical properties. 

In the current work, we formally assume that the Reynolds averaging operator is a 
linear operator, conventionally denoted by an overbar, with the following properties: 

(1) Linearity: for two variables hi and h 2 , and two numbers ai and a 2 , 


(3.1) 


aihi — aih\ 


Additionally, a = a if a is a number. 

(2) Commutativity with respect to differentiation and integration: 


(3.2) dh = dh, j h = j h. 

(3) The Reynolds property: 

(3.3) hih2 = hi ii2- 


We also formally assume that the Reynolds-averaged quantities p, u, p, as well as their 
Reynolds-averaged products (such as for example) are sufficiently smooth at small 
scales to be computationally tractable. However, due to inherent nonlinearity in the ve¬ 
locity variable u, the transport equations for these Reynolds-averaged variables involve 
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new unresolved variables which are the Reynolds averages of products of turbulent fluc¬ 
tuations between the fine-grained and averaged variables. For example, the averaged 
momentum equation involves the Reynolds-averaged outer product of the turbulent ve¬ 
locity fluctuations, called the "Reynolds stress", for which a closure is usually sought 
in the form of the Boussinesq approximation |^. In this work we present a different 
approach, which produces the equations for the coarse-grained variables with a better 
closure potential. 

Here we formally apply the Reynolds operator, described above, to the Boltzmarm 
equation in ( 2A) (rather than the Navier-Stokes equations), as was done previously by 
Girimaji | [T^ . The main novelty of our work is that, unlike p^ , we further rebuild 
the hierarchy of the transport equations like the one in ( |2.7[ > for the moments of the 
Reynolds-averaged Boltzmarm distribution. The resulting coarse-grained moment vari¬ 
ables are different from the usual Reynolds averages of conventional moments in ( 2.5| >, 
and we provide suitable formulas to relate the new coarse-grained moment variables to 
the Reynolds averages of the conventional moments. The key advantage of our approach 
is that the transport equations for the new coarse-grained moments are identical to the 
conventional moment transport equations in ( |2.7| > (collision terms being an exception), 
and, therefore, naturally lend themselves to the variety of previously developed closures 
for (2.7). The only unspecified closure parameters originate from the nonlinear collision 
term of the Boltzmarm equation, which is not explicitly tractable by a linear averaging 
operator. 

We define the new coarse-grained variables p*, u*, p*, S*, q*, Q*, R* as follows: 


(3.4a) 

(3.4b) 

(3.4c) 

(3.4d) 

(3.4e) 

(3.4f) 


p* = (l)j, coarse-grained density. 


p u = {V)j, 

1 


coarse-grained momentum. 


p* =-(||z; — M*II )j, coarse-grained pressure, 

S* = {{v — u*) ^ {v — u*))j — p*I, coarse-grained stress. 


1 


— 2^11^ “ ^*lr(^ “ coarse-grained heat flux, 

Q* = {{v — u*) ® {v — u*) ® {v — u*))j, coarse-grained skewness. 


(3.4g) R* = -(||u — M*||^(u — M*) 0 (u — u*))j, coarse-grained contracted flatness. 
Then, the following straightforward result follows. 

Proposition 2 (Coarse-grained transport equations). The new coarse-grained variables are 
related to the Reynolds averages of the fine-grained variables as follows: 

(3.5a) p* = p, 


(3.5b) 


- 

p u = pu, 


u' = u — u*, 
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(3.5c) 

(3.5d) 

(3.5e) 


p* = p + giollw'p/ 
1 - 


S* = S + pu^ ®u' — -pWu^W^I, 

o 

^ 5—, 1- 


q* = q + Su' + '^pu' + ^p\\u'\\^u', 


TT 


(3.5f) Q* = Q + (pi + S) ^u'+ (pi + S) ® u' +{pI + S)^u' +pu'^u'®u', 


(3.5g) R* = R + Qu' + [Su' + q) ®u' + u' ® (Su' + q) + 


+ -pu' 0 m' + - \\u'\\^(pl + S) + -p\\u'W^u' ^ u'. 

Above, and A^^ denote two possible transpositions of a 3-rank tensor A. The transport 
equations for the coarse-grained variables in ( |331 > are of the same form as the moment transport 
equations in namely 


(3.6a) 

(3.6b) 

(3.6c) 

(3.6d) 


d{p*u* 

dt 


^ + div{p*u*) = 0, 


+ div{p*u* 0 M* + p*I + S*) = 0, 


dS* 

dt 


+ div(p*M*) + ^ [p* divu* + S* : (V 0 u*) + div^*] = 0, 


+ (u* ■ V)S* + div(M*)S* + 2 [S*(V 0 u*)] 


+ divQ* - -(div^*)I + 2p* [V 0 u*] = Cg, 


(3.6e) f + div(,- ® „•) + (,• . V)„- - 1 (Ifl + S-) div(p-I + S-) + 

+ Q* : (V0M*) + divR* = 4 

where the coarse-grained collision terms are given by the Reynolds averages 
(3.7) cj = + Csu'. 

Proof The transport part of ( |3.6| > is identical to that of ( |2.7| > since the transport parts of 
the underlying Boltzmann's equations for f and / are identical. To derive the formulas 
in ( |3.5| >, first observe the following relafions: 

(3.8a) p* = {l)j = Wf = ^ 

p*u* = {v)j= {v)f = pu, 


(3.8b) 
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(3.8c) p* + ^p*\\u*f = ^{\\vf)j= l(||x;||2)^ = lp||„||2^ 

(3.8d) S* + p*I + p*u* = {v® v)j = {v ® v) f = S + pi + pu ® u, 

(3.8e) + S*M* + = '^{\\vfv)j = 

1 - _ _ 5 _ 1 - 

= 2^11^11^^)/ = ^ + Sm+ 2?" + 2^11^11 V 

(3.8f) Q* + (p*I + S*) ®u* + {{p*I + S*) ® u*f + {{p*I + S*) (g) u*f^+ 

+ p*u* ®u* ®u* = {v®v® v)j = {v®v®v)f = Q+ 

+ {pi + S) 0 M + {{pi + S) 0 uY + {{pi + S) 0 uY^ + puWuWu, 


(3.8g) R* + Q*u* + {S*u* + q*) ®u*+u*® {S*u* + q*) + aP*"* ® «* + 


+ ^\\u*Y{p*I + S*) + ^p*\\u*Yu* ®U* = ^{\\vYv®v)j= ^(||z7||2t7 0Z7)^ = 

1 - 


1.2, - - 1- 


= R + Qu + {Su + q)®u + u® {Su + q) + -pu ®u + -\\uY{pI + S) + 

1 - 


+ -| 0 ||m|| 2 m 0 u. 


Rearranging the above relations and noting that pu' = 0, one obtains the definitions for 
the coarse-grained variables in ( 3.5| >. Lastly, let us derive the formulas for the coarse¬ 
grained collision terms. The Reynolds-averaged moment transport equation in ( 2.3| > is 
given by 


(3.9) 


— {b)j + divYbv)j= (b)^. 


In particular, for the stress matrix it is given by 


(3.10) 




v)j+divx^ 


V ®V ' 


'Yl = Csr 


where the collision term is unchanged, since the combinations of equations, added to the 
stress equation, have no collision terms. Now, subtracting the same combinations of the 
transport equations for the coarse-grained variables (which also do not have collision 
terms) in reverse order, we obtain the coarse-grained stress equation with the same 
averaged collision term. For the heat flux equation the situation is slightly different, as 
it is given by 

13 1 

(3.11) -^^{\\vYv)j+-dwx{\\vYv®v)j = + Csu, 


as the derivation of the moment equation above requires adding the stress equation, 
multiplied by u, to the heat flux equation. Now, the derivation of the coarse-grained 
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heat flux equation from the equation above requires subtracting the coarse-grained stress 
equation, accordingly multiplied by the coarse-grained velocity u*. This results in the 
coarse-grained collision terms of the form (3.7|>. □ 


Unfortunately, there appears to be little we can do about the Reynolds-averaged colli¬ 
sion terms in the context of turbulence transport. While the transport part of the Boltz- 
maim equation is linear and thus is "transparent" to the Reynolds averaging, the same 
is not true for the collision terms. Without more information about the properties of the 
statistical ensemble (and, therefore, statistical physics of the fluid under consideration), 
we can only conclude that the coarse-grained collision terms C 5 and have a general 
dissipative effect towards the global maximum entropy state under the total mass, mo¬ 
mentum and energy constraints. Below in Section]^ we suggest crude parameterizations 
of Cg and c* based on dimensional arguments. 


4. The transport of turbulent energy under the anelastic turbulence 

APPROXIMATION 


In the previous section we derived a hierarchy of the coarse-grained equations ( |3.6| >, 
which does not seem to require anything other than a suitable closure to become com- 
putationaly tractable. However, observe that the new coarse-grained variables in ( 3.5| > 
are not exactly what one would naturally need to model. For example, a meteorologist 
would likely need to model, among other variables, the Reynolds average of the pres¬ 
sure p (as the closest computable alternative to the pressure p itself), yet observe that 


the coarse-grained variable p* is not p, but rather p + pWu'W^/3 (see ( |3.5c[ )), where the 
additional term is not zero! 

Therefore, we need to provide a means to recover the Reynolds averages of the fine¬ 
grained variables from the computed coarse-grained variables. In order to do this, we 
make an additional simplifying approximation about the nature of the small-scale tur¬ 
bulent fluctuations, which absolves us of the need to address the Reynolds averages of 
the density and velocity variables, leaving only the pressure, stress and heat flux to work 
with. Above, in the derivation of the equations for the coarse-grained variables, we as¬ 
sumed the most general form of turbulence, were the turbulent fluctuations were present 
in all fluid variables, and thus the Reynolds averages of fluid variables were generally 
never identity operators in the presence of turbulence. However, in many applications it 
is known that the density p is weakly affected by the small-scale turbulent fluctuations, 
and varies only on the coarse-grained "resolved" scale, while the higher order moments, 
such as the velocity u, pressure p, etc are still affected by the turbulence. Thus, here 
we assume the identity p = p*, and, therefore, p can be factored out of the Reynolds 
averages. As a result, the turbulent velocity component u' = u — u* becomes anelastic: 

(4.1) div(|OM') = 0. 


Indeed, observe that subtracting the coarse-grained density equation ( |3.6a| > from the 
fine-grained density equation (2.7a|> results in 


(4.2) 
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Setting p* = p above yields ( |4.1[ ). 

The anelastic turbulence approximation does not change the form of the coarse-grained 
transport equations in ( |3.6| >, but instead simplifies the definitions of the coarse-grained 
variables in ( 3.5| > to some extent. For the Reynolds averages of the anelastic turbulent 
fluctuations 


(4.3) 

we obtain directly 
(4.4a) 


u' = u — u*, 


V' = V- V*r 


S' = S- S*, 


q' = q- q*, 


u' = 0, 


(4.4b) 


v' = -\p\WV> 


(4.4c) 


1 , 


s' = —p [u’ ®u’ — -llw'lpl j , 


(4.4d) 


1 


q' = — [ S'u' + -p'u' + -p\\u'\\^u' ) . 


Obviously, computing p'. S' and q' in addition to the coarse-grained variables in ( |3.5| > 
allows to coimect back to the Reynolds averages of the original variables. Here we are 
going to start with p', which requires its own separate equation. 

To derive the transport equation for the Reynolds-averaged turbulent pressure p', we 
subtract the coarse-grained pressure equation ( 3.6c| > from the fine-grained pressure equa¬ 
tion (2.7c|> and obtain 


(4.5) ^ div {p*u' + p'u* + p'u') + ^ (p^divw' p'divw* p'divw') 

^ (S : (V O m) — S* : (V O m*)) ^div^' = 0. 
Applying the Reynolds average to the equation above results in 


dp' 


2 ^ 


(4.6) + div yp'u*J + -p'divu* + div i^p'u'j + -p'divu'+ 


+ - yS : (V 0 u) — S* : {V u*)j + -div^' = 0 . 


2 ^ 


We now observe that 


(4.7a) p'divw' = —pdivu' = —(div(pu') — u' ■ Vp) = —p'u' ■ 

p p p 


(4.7b) 


S : (V 0 m) - S* : (V ® M*) = S' : (V ® M*) + S' : (V 0 u'), 
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which further results in 


ap' 


2 ^ 


(4.8) + div ( p'm ) + -p'divM + div ( p'w') — - • p'u'+ 


2V£.;q-7. 


3 p 

+ ^ (^ : (V ® M*) + S' : (V0mO + div^) = 0. 


The turbulent pressure p' is not a convenient measure of turbulence from a physicist's 
perspective, since, according to ( |4.4b| >, it caimot be positive. Thus, it is more convenient 
to replace it with a noimegative quantity, and here we choose the kinetic energy of 
turbulent velocity fluctuations |25 26| for this purpose: 


(4.9) 


Et = -\p' = \p\\u'V. 


For the turbulent energy in ( |4.9| >, the transport equation becomes 


(4.10) + div(EiM*) + ^Etdivu* — ^div(p'M') + • p'u' — 

C/t vj ^ 


— s' : {V 0 u*) — s' : {V u') — divq^ = 0. 


The term : (V (8) u') in the equation above can be rearranged as follows. We assume 
that the microscale viscosity p can be factored out of the Reynolds average (that is, 
replacing p with its own Reynolds average is an acceptable simplification). Then, one 
writes 


(4.11) S' : (V <S) u') = S : (V <S) u') = —2p [V (8) w] : (V (8 u') = 


= -2p [V 8 w'] : (V 8 u') = -p(divM')2 — p(V 8 u') : (V 8 u'Y — p|| V 8 u'\\ , 


where we replaced the fine-grained stress S with its Navier-Stokes approximation ( 2.27a| >. 
For the first term in the right-hand side of (4.111 we use the anelastic turbulence approx¬ 
imation (14.11) and obtain 


(4.12) {divu'Y = ^(|0divM')2 = ^(div(pM') — Vp • u'Y = ® ® ^ 

= -(S' + fl) : i(Vp®Vp) = ?l!^E,-l(Vp®Vp) 
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where we used ( 4.4c| > and (4.9) in the last equality. For the second term in the right-hand 
side of (4.11|) we obtain 


(4.13) (V 0 u') : (V 0 u'Y = (V 0 V) : {u' 0 u') + 2div (^u' 0 + 




+ ^(V |0 0 V| 0 ) : {u' 0 m') = (V 0 V) : ^Etl - S' + 2div ;^EtVp - + 


3p 


d'^P 


1f7 


3p 


+ -^(Vp 0 V| 0 ) : —Efl - -S = -A — + -div Ef^ + - 


Et 


3p^ 

Vp\ 2||V|0| 


S' 


d^P 


- (V 0 V) : - - 2div S'^ {\/p®Vp): S' = - 


3 p3 
^_2AEt , 2^Ap 


-Ef- 


+ T^Et ^— 




3 p 3 p 

IZ^_1(V0V) :^-^((V0 V)p) :¥+^(Vp(8)Vp) 


Now, we denote the last term in the right-hand side of (4.11) as the turbulent energy 
dissipation rate 125 26 381, 


(4.14) 


pet = ^11V 0 u'W , 


as it caimot be expressed in terms of the Reynolds averages of the turbulent fluctuations 
that we already defined. This leads to 


(4.15) S': (V0 mO = -l^AEt + l^ 


3p 


+ ^(V0V) :S'+^((V0V)^i) : S'-^(Vp0Vp) : S'- pe^ 


3p \3 p^ 

p 


P J 


5p 


and, consequently. 


(4.16) ^ + div(EfM*) + ^ l^divM* + Ef- 

- (^V0M* + ^(V0V) + ^((V0 V)^) - ^(Vp0Vp)^ :W- 

— div^' — -div • p’u' -E pet = 0. 

Observe that the equation for the turbulent energy above depends on the Reynolds- 
averaged turbulent stress S', heat flux q', the quantity p'u', as well as the turbulent 
energy dissipation rate pet- 


4.1. Isotropic turbulence assumption and its limitations. As a special case of the tur¬ 
bulent flow, here we consider a hypothetical situation where the turbulent velocity fluc¬ 
tuations u' at each spatial point are statistically decorrelated and isotropic. This implies 
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the following conditions for statistical averages of turbulent quantities: 


(4.17) 


u' ®u' = - ||w^||^I, 

O 


s' = 0, Tp'u’ = q' = Q. 


The first equality above is due to the fact that the cross-correlations between different 
components of the turbulent velocity are zeros, while self-correlations are equal. The 
second equality follows immediately from the first one by using ( 4.4c| >. The third equality 
is due to the fact that, because of isotropy, for each instance of u', p' and q', there must 
be another with —u', p' and —q', and with the same statistical weight. This is, of course, 
an idealized assumption, which is unlikely to hold exactly in practical situations. In fact, 
it is somewhat of a turbulent analog of the Maxwellian equilibrium condition in ( 2.12| > 
for the fluid particle velocities, assumed without rigorous justification. 

Under the isotropic assumption in ( 4.17| >, the equation for the turbulent energy trans¬ 
port in (4.16|> becomes 


(4.18) ^ div{Etu* 


+ ^ ^divM* 


M . Ap 5 
“ 3 ''- 


Vp| 


Et pst — 0 . 


There are two problems with this equation. First, the diffusion term above is clearly 
ill-posed, which means that the turbulent energy will be amplified by this term at small 
scales, rather than damped. The second problem with the equation above is that it has 
no means of producing turbulent energy from zero initial condition. At the same time, 
it is very well known that the turbulence appears even if there was none to begin with, 
from the interactions of the large-scale motions with the small-scale viscous dissipation. 
Thus, the only way to create such interaction is to couple the higher-order Reynolds 
averages in ( |4.16| > to the coarse-grained variables, rather than assuming isotropy. In a 
broad sense, the turbulence is produced by the anisotropy of the large-scale flow. 


4.2. Coupling the turbulent Reynolds averages to the coarse-grained variables. Here 
we assume that on the fine-grained scale the fluid satisfies the Navier-Stokes approxi¬ 
mations in ( 2.27| >. We also assume that replacing the viscosity p and heat conductivity 
K in ( 2.27| > with their corresponding Reynolds averages is an acceptable approximation. 
Then, the Reynolds averages of the turbulent fluctuations S' and q^ can be expressed via 
the coarse-grained stress S* and heat flux q*, and the Reynolds-averaged Navier-Stokes 
approximations for the fine-grained stress and heat flux as follows: 


(4.19a) 


S' = S-S* = -S*-2p[V' 


>u*] 


TS' 


cj'= =-cj*- kV { = -q* -kV (^j +^kV 


(4.19b) 


P 
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Substituting the above expressions into the transport equation for the turbulent energy 
yields 


(4.20) 


^ + div{Etu*) + ^ (divu* + yA + 


3 V 


Et — -div ( kV ( — 


+ 


+ (v ®u* + ®V)p) - ^{Vp ®Vp)^ : + 


+ 


+ div^* + div ( kV ( ^ J j — ^div(p'M') + • p'u' + pet = 0. 


Vp 


The unknown term p'u' is still present, and needs either its own separate transport equa¬ 
tion (which will, of course, involve yet higher-order moment combinations of turbulent 
quantities), or a suitable closure. For simplicity, here we propose a closure under the 
assumption of Bernoulli's principle. For that, first observe the following relation: 


(4.21) S'u' = Su' = —2p [V (8) u\ j^u' = —2p [V (8 j^u' = 


= -p 8 w' (V 8 u'Y — -(divM')7j u', 

where, as before, the Navier-Stokes approximation is used for the Reynolds average of 
the fine-grained stress S. For the separate terms above we have, with the help of (4.1), 


(4.22a) 


(V8«0«' = -yV(^) =V(^), 


(4.22b) 


(V8wOV = 


1 

P 


div(pM' 8 u') 


1—2 
—divS' + - V 
p 3 



2p Vp 

3^^ p2 ' 


(4.22c) 


(divM')M' = —(«' 8 u')^-^ = 


p p2 3 


which, when assembled together, lead to the expression of S'u' in terms of the Reynolds 
average of the turbulent stress fluctuation S' and the turbulent kinetic energy Et: 


(4.23) 


F. 


2p 


F7 5, 


S'u' = ( ^div + ^Vp ) S' - -pV ( ^ 


^p 3p2 '^J~ 3 

Combining (4.23) with (4.4d) and ( |4.19| ), we obtain 


(f) 


10 Vp 


(4.24) 
q* + KV 






div -|- 



(s*+2p[V8w*]^s) - 



lOp 

9p2 


EtVp + -p'u' + -p\\u' 


\^u' 
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or, after rearranging the terms. 


1 


(4.25) -p'w' + -| 0 ||m'Pm' = q + kV I — + 


P 


5}i — 2 k 




V — 
\P 

2p 


+ 


lOp 

9p2 


EiVp+ 


Observe that the Reynolds averages of turbulent fluctuations in the left-hand side above 
are expressed entirely in terms of the coarse-grained variables and the turbulent energy 
Et in the the right-hand side. However, in order to close the turbulent energy transport 
equation in terms of the coarse-grained variables and the turbulent energy itself, one has 
to surmise a plausible relation between the Reynolds averages p'w' and There is 

certainly more than one way to do that (as an example, one could postulate = 0, 

for a simplest closure), and we here choose what we perceive as the most physically 
realistic. Namely, we recall Bernoulli's principle, which states that the change in the 
local flow velocity tends to affect the pressure in the flow in a manner as to preserve the 
total energy. If we treat the turbulent velocity fluctuation u' as the primary cause for the 
turbulent pressure fluctuation p' at the same location, Bernoulli's principle then leads to 
the relation 


(4.26) 


/ 1 II /||2 

P =-3Pll«ll . 


Of course, we understand that in reality the relation in ( |4.26| > caimot hold at each spatial 
point exactly, as Bernoulli's principle is merely an approximation. However, we note 
that the exact relation ( 4.4b| > is precisely the Reynolds-averaged relation ( 4.26| >. So, we 
surmise that it is plausible to incorporate an approximate relation in ( 4.26| > into the 
relation between p'u' and as follows: 


(4.27) 


1 


o 


or 


-p'w' -t -p||m'Pm' = p'm'. 


In particular, as a result of Bernoulli's assumption, the Reynolds average of the turbulent 
heat flux becomes the Reynolds average of the product of the turbulent velocity with the 
turbulent pressure matrix, taken with the opposite sign: 


(4.28) 


= -p'u', P' = p'l + S'. 


This results in the following expression for p'u' in terms of the coarse-grained variables 
and the turbulent energy Ef. 



+ 


5p — 2 k 
3 


V 



+ 


lOp 

9p2 


El Vp-|- 


+ 



div -|- 



(s*+2p[V®M*]^s) . 


(4.29) p'u' = q* + kV 
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The turbulent energy equation then becomes closed with respect to the coarse-grained 
variables and the turbulent energy Et itself: 


(4.30) 


^ + div(EtM*) + ^ (divu* + ^ Et - ^div (kV 


3 V 




+ (- ?VV f f + (1^ 4V kV f ^ 1 1 + 


+ ( V®M* + ^(V®V) + ^((V®V)p) -^(V^®Vp)) : [S* + 2}i[V0u*]^^) + 


n 


5}i 




^div+4vp 

P 


S* + 2}i[V®u*] 


TS 


+ pct — 0. 


Rearranging the terms, we further obtain 


(4.31) 
dEt 


^ + div { Et [u* + 


15}i — 4 k V(2k —15^) 

i2 ‘ 


3p^ 


6p 




+ 3 "^ 




+ (K-4^)^ + V(K-5^)-4 + divf^^4r4^1 1 = E; 




4p 




6p 


pet, 


where is the coarse-grained forcing: 


(4.32) Ef = ( ^ 


^* + kV(4 )- 


(^®u* + ^{V®V) + ^{{V®V)p)-^{Vp®Vp)^ : (s* + 2 ^[V 0 m*]^ 5 ) + 

'S* + 2^ [V (8) w*] ^5 




^^div + 

P 3p2 


Here observe the following properties: 

• For the Prandtl number of the ideal gas, Pr = 2/3, and the assumption that the 
spatial derivatives of p and k can be neglected, the advection term in ( 4.31| > is 
exactly the same as it is for the usual pressure variable. 

• The Laplace diffusion term in ( 4.31| > is well posed as long as Pr > 5/11, which 
includes the ideal gas. 

• The term with || V| 0 |p in (4.31) is linear damping as long as Pr > 1/3, which again 
includes ideal gas. 

For a strong turbulence (that is, large coarse-grained stress S* and heat flux q*) and small 
microscale viscosify p and heat conductivity k, one likely can drop the terms scaled by 
the latter in ( |4.31[ ) and ( |4.32| >, except the one for the diffusion, as it could potentially be 
the highest-order differential operator, given the lack of information on pe^. The resulting 





















COARSE-GRAINED TRANSPORT VIA MOMENTS OE AVERAGED BOLTZMANN EQUATION 


23 


simplified equation is given by 
(4.33) 

^ + div{Etu*) + ^Etdivu* + = QV - - S* : (V (8) «*) - pet- 


Observe that this simplified transport equation tor the turbulent energy is the same as 
the one for the Euler approximation of the pressure equation, but with the additional 
forcing via the coarse-grained variables S* and q*, small scale dissipation pst, and a 
weak diffusion. 

At this point, it becomes clear how the turbulence is produced from an initial condition 
which belongs to the fully resolved coarse-grained scale, with zero stress and heat flux. 
Schematically, this process can be illustrated as 


(4.34) 


2p*[V®u*] 


TS 


A q^ 


Et, 


where the are two stages: 

(1) As the stress and heat flux are zero initially, the coarse-grained strain rate term 

2p* [V and the coarse-grained temperature gradient term {5/2)p*'V{p* /p) 

act as external forcing in the equations for the coarse-grained stress ( 3.6d| > and 
heat flux ( 3.6e| >, respectively, causing deviation from the Maxwellian equilibrium 
on the coarse-grained scale. 

(2) In turn, the nonzero coarse-grained stress S* and heat flux act q* act as external 


forcing in the turbulent energy equation (4.31), causing the nonzero turbulent 
energy to appear. 

In particular, this means that the coarse-grained collision terms Cg and provide 
weaker damping rate than their tine-grained counterparts even when there is no tur¬ 
bulence developed yet. 


5. Collision terms and moment closures for the coarse-grained transport 

EQUATIONS 


Thus far, we developed the following ingredients of the coarse-grained transport 
framework tor a turbulent flow: 






A hierarchy of the new coarse-grained variables in (3.5) and the corresponding 
hierarchy of the transport equations in ( 3.6| ). Although the expressions for the 
coarse-grained variables in ( 3.5| ) and the corresponding equations ( 3.6| ) only in¬ 
clude the moments up to the heat flux, clearly the full set of moments is infinite, 
and requires either infinitely many transport equations, or an approximate clo¬ 
sure. The derivation was done under the general assumption that, microscopi¬ 
cally, the flow is described by the Boltzmaim equation for a monatomic ideal gas 

(|g. 


The relations for the corresponding Reynolds-averaged turbulent fluctuations. 
Those are given in the form of the turbulent energy (or, equivalently, turbulent 
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pressure) transport equation in ( 4.31| >, and the relations for the turbulent stress 
and heat flux in (4.19). This was done under the following assumptions: 

(1) Validity of the Navier-Stokes approximations for the fine-grained stress and 
heat flux (12.271); 


(2) Invariance of the density p under the Reynolds averaging and the resulting 
anelastic turbulence approximation ( |4.1[ ); 

(3) Approximation of the microscale viscosity p and heat conductivity k by their 
own Reynolds averages; 

(4) Bernoulli's principle for the turbulent velocity fluctuations in ( 4.27[ ). 

In order to make the resulting system of equations suitable for practical computations 
and modeling, two further parameterizations need to be developed: 

(1) A closure of the hierarchy of the coarse-grained equations in ( 3.6| ). Namely, we 
need to relate the coarse-grained skewness moment Q* and the contracted fourth- 
order moment R* to the coarse-grained density, pressure, stress and heat flux. 

(2) Suitable parameterizations of the coarse-grained collision operators Cg and c*, 
and the turbulent energy dissipation rate pet- 

Below we elaborate on the first subject to the extent it allows to avoid the second subject. 
The reason for this is that the collision operators and the turbulent energy dissipation 
rate are the Reynolds averages of nonlinear (with respect to the transported variables) 
quantities, and thus should likely be modeled depending on a particular application, 
with more detailed assumptions on the flow properties put into place. Nonetheless, in 
Sectionj^we suggest some crude parameterizations for Cg, c* and pet. 

5.1. Grad closure for the coarse-grained transport equations. Generally, there is little 
that can be concluded about the local entropy state Si[f] of /. As Girimaji pointed out 
in 1121, the problem here is that / is the statistical ensemble average of many realizations 
of /Twhere each realization of / generally has its own density p, velocity u, and pressure 
p at a given spatial point x. Therefore, even if we assume that each / in the statistical 
ensemble is the corresponding Maxwellian state ( 2.12| ) with its own p, u and p, the 
averaged sum of these states does not have to be near a Maxwellian state. However, if f 
is far from a Maxwellian, then both the Euler and the Navier-Stokes approximations for 
the moments of such distribution become rather questionable. 

In contrast to the local entropy state, it turns out that the global entropy state of the 
Reynolds-averaged Boltzmann equation behaves similarly to the Reynolds average of the 
global entropy. In order to better argue our point here, we first note that the Reynolds 
operator satisfies Jensen's inequality. 

Proposition 3 (Jensen's inequality). The Reynolds operator satisfies Jensen's inequality: for a 
convex function <p and a function of the statistical ensemble h, 


(5.1) (p{h) > (p(h). 

Proof Follows from the fact that the Reynolds operator is a statistical average. 


□ 


This results in the following estimate for the entropy of the Reynolds-averaged solu¬ 
tion of the Boltzmaim equation (|2.1 ): 
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Proposition 4 (Entropy of the Reynolds-averaged Boltzmann equation). Let f be the 
Reynolds-averaged distribution f of the Boltzmann equation in which apparently satis¬ 
fies 


(5.2) -^ + vV,f = C(J). 

Then, both the local and global entropies of f are bounded from below by the Reynolds averages 
of the local and global entropies of f, respectively: 

(5.3) S,(7] > W], Ss[71 > Wl- 

Proof Observe that (p{f) = / In / is convex. Therefore, Jensen's inequality yields 

(5.4) Si[f]{t,x) = - Jflnfdv > - J /In/du = S//](f,x). 

The corresponding inequality for the global entropy is obtained by integrating in x. □ 


The above result is not of much use for estimating the local behavior of / at a spatial 
point x; indeed, while the lower bound estimate is valid, the upper bound for Sff] can 
be different from the upper bound for Sff], due to the fact that the local density, mo¬ 
mentum and pressure are not preserved between the members of the statistical ensemble 
(if they were, there would likely be little need for the Reynolds averaging to begin with). 
However, it is quite reasonable to assume that the members of the statistical ensemble 
share the total mass, momentum and energy constraints, because they have to be dif¬ 
ferent realizations of the same large-scale flow scenario, rather than being completely 
unrelated. Therefore, the statistical ensemble members also share the maximum global 
entropy state under these constraints, which means that the maximum global entropy 
state is invariant under the Reynolds averaging. As a result, one can show the following 
inequality for the global relative entropy, 

(5.5) Hg[f] = -Sg[f] - Jfhxpl^dvdx < -Sg[f] - J flnpl^dv dx = Hg[f]. 


The latter means that the global relative entropy of the Reynolds-averaged distribution 
/ decays at the same rate as the Reynolds average of the ensemble global entropy states 
(that is, 0(f““), as shown by Desvillettes and Villani |10|). 

This situation is similar to what usually happens in the dynamics of rarefied gases 
p 14 T5|2T|22| where there is little time scale separation between the growth rates of the 
local and global entropies, and thus the higher-order moments do not rapidly converge 
to their local equilibrium values, while still maintaining global convergence on a slower 
time scale. Thus, the Grad closure appears to be a suitable option tor closing 

the transport equations in ( |3.6| > for the coarse-grained variables in ( |3.5|>. Applying the 
Grad closure to the coarse-grained moment transport equations in ( |3.6| > is identical to 
that of the usual moment transport equations in ( |2.7| >, as the Grad closure imposes no 
restrictions on the structure of the collision terms, whether Reynolds-averaged or not. 
The resulting system of equations for the coarse-grained variables is the same as the one 
























26 


RAFAIL V. ABRAMOV 


in ( 2.32| >, except that the collision terms are different; 

^ + div(pW) =0, 


(5.6a) 

(5.6b) 

(5.6c) 

(5.6d) 

(5.6e) 


d{p*u* 

dt 


+ div{p*u* O M* + p*I + S*) = 0, 


+ div(p*M*) + ^ [p* divu* + S* : (V (8) u*) + divq* 


= 0 , 


ds* 

dt 


dt 


+ (m* • V)S* + div(M*)S* + 2 [S*(V (8) u*)] + - [V ® q*] 

+ 2p*[V®M*]^s = q. 


+ -(div(g^ 


u*) + {q*-VX) + -{iV~ 


>u*)q* - {u* ■V)q*)- 




div{p*I + S*) + div 


5p*%_l_7p* 


s* = c! 


2p* Ip’* 

The following arguments can be made in favor of the Grad closure in ( 5.6| >: 

• The ability to set the prescribed stress and heat flux. It caimot be assumed that 
the coarse-grained stress S* and heat flux q* are zero (or nearly zero), as there is 
not enough information about rapid Reynolds-averaged local entropy growth (or, 
equivalently, decay of Hj [f]) to substantiate that. Thus, a suitable moment closure 
must account for that. At the same time, the moments of such closure must be ex¬ 
plicitly computable in terms of elementary functions, for the transport equations 
to have an explicit form. Also, one must remember that, due to Proposition the 
Reynolds average / has a global tendency to dissipate towards a normal distribu¬ 
tion on a long time scale, and thus the closure should not be "too far away" from 
this state. The Grad closure meets all of these conditions: it is built around ( 2.12| >, 
its moments are explicitly computable, and it also possesses the prescribed stress 
and heat flux, in addition to prescribed density, velocity and pressure. 

• Few restrictions for the parameterization of the coarse-grained collision terms. 
As the coarse-grained collision terms Cg and are nonlinear and require sep¬ 
arate (quite possibly empirical or semi-empirical) treatment, any need in prior 
assumptions on the structure of the collision terms (as in the Navier-Stokes ap¬ 
proximation, for example) is undesirable. In the Grad closure, there are no a 
priori assumptions on the form of collision terms, as they enter the equations as 
unspecified parameters. 

• Backward consistency with a weak/vanishing turbulence scenario. If the turbu¬ 
lence is weak or vanishing, then the coarse-grained variables in ( |3.5| > become the 
corresponding fine-grained variables in ( |2.5| >. Then, the Grad equations for the 
coarse-grained variables naturally transition into the ordinary Grad equations in 
( |2.32| >, with appropriate collision terms. 
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• Generalization of Millionschikov's hypothesis. In 1941, Millionschikov 
asserted that the relations between the second and fourth turbulent statistical 
moments are related in approximately the same way as the corresponding mo¬ 
ments of the normal distribution. Later in 1948, Heisenberg independently 
proposed the same hypothesis. Since then, some evidence was accumulated in 
favor of fhis hypothesis (see |2^|26| and references therein). The Grad distribu¬ 
tion generalizes Millionschikov's hypothesis, providing the explicit dependence 
of higher-order turbulent moments on the stress and heat flux in addition to the 
density, velocity and pressure. 

However, there are also the following drawbacks: 

• Increased computational complexity. Observe that the usual Navier-Stokes equa¬ 
tions incorporate only five prognostic variables: the density, three velocity com¬ 
ponents, and pressure. For the Grad closure, the number of prognostic variables 
increases to 13 (plus one more for the turbulent energy Et). 

• Somewhat questionable numerical stability. There are at least two types of pos¬ 
sible numerical instabilities which can manifest in the Grad equations. The first 
one arises when the contribution of the collision terms for the stress and heat flux 
is very large, so that these variables are strongly damped, and thus the numeri¬ 
cal stiffness (thaf is, oscillation of the numerical time-discretization polynomial) 
arises for high-accuracy integration schemes unless the time step is very small, or 
a specially tailored scheme, such as an implicit BDF, is used. The second numeri¬ 
cal instability will arise if the collision terms Cg and do not include diffusion, 
in which case the numerical solution will develop shock waves |14| which will 
transfer the energy to high Fourier wavenumbers and cause the Gibbs oscilla¬ 
tions. Below in Section]^ we suggest an option to mitigate this problem. 

5.2. Coarse-grained viscosity and heat conductivity approximations. Observe that as 
the turbulence vanishes, the coarse-grained variables in ( 3.5| > naturally become the fine¬ 
grained variables in ( 2.5| >, while the coarse-grained transport equations in ( 3.6| > become 
the fine-grained transport equations in (|2.7|>. The latter suggests that, under the condi¬ 


tions of sufficiently weak turbulence, the coarse-grained collision terms Cg and can be 
represented as linear damping via the coarse-grained viscosity ji* and heat conductivity 
K* as follows: 

(5.7a) 


Cs 


——S* 


(5.7b) 


/-r ' _ 

~ - Z" 


2k*^ ' 


where and k*, of course, remain unspecified quantities and must be somehow pa¬ 
rameterized, possibly semi-empirically. However, this situation is much simpler than 
what it was before from the modeling perspective, as now two scalar quantities need to 
be estimated, rather than two general collision operators. Below in Section we sug¬ 
gest a crude approximation for both ji* and k*, which could serve as a starting point in 
practical modeling, at least until the properties of Cg and c* become better studied. 
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The main drawback of the linear damping parameterization is the following. In prac¬ 
tice, the microscale viscosity and heat conductivity are very small, while the turbulent 
fluctuations are large. In a typical turbulent flow situation, S' ^ S and (j' S> (j, such that 
the coarse-grained stress S* and heat flux q* consist largely of the Reynolds-averaged 
turbulent fluctuations. 


(5.8) 


S* ^ -S', 




-q'. 


In this setting, there is simply not enough information to determine whether pS and 
^ (which are very small quantities) are collinear, respectively, to p*S* and p*q* (which 
are large quantities), as the latter are determined largely by the statistical physics of the 
turbulent flow, rather than microscale collisions. However, this collinearity is a necessary 
condition to justify scalar quantities p* and k*. 

5.3. Coarse-grained Navier-Stokes equations. If, in addition to ( |5.7| >, the coarse-grained 
viscosity p* and heat conductivity k* are sufficiently small, it enables the Navier-Stokes 
approximation of the form ( |2.25[ )-( [Z^ , that is, 

(5.9a) S* =-2p* [V ® u*] 

(5.9b) q* = -K*VipVp), 

so that the coarse-grained Grad equations in ( |5.6| > simplify to the coarse-grained Navier- 
Stokes equations: 


TS' 


(5.10a) 

(5.10b) 


^ + div(pW) =0, 


d{p*u* 

dt 


+ div {p*u* ^u*) + Vp* =2 div (^p* [V ( 8 ) u*] 


(5.10c) ^ + div(p*M*) + divM* = ^p* I [V G) u*] ^ J + ^div(K* V(pVp*)), 

where, as in ( |2.28| >, one could consider dropping the term with p* in the coarse-grained 
pressure equation, since it acts as unbalanced forcing which tends to increase p*. In this 
situation, the coarse-grained forcing in the turbulent energy transport equation ( 4.31| > 
becomes 


(5.11) F,* = I N - ^ 


K — K 


'>^(7 


-2{p-p*) (v + ^{S/ ®V) + ^V)p) - ^(S/p^S/p)) : 

\ p p 3p J 


TS 


+ 


yp 


+ ( 3V - 2yF ) . ( ( ^div + ) ({p - p*) [V » »•] 


2p 


TS 


Assuming that the coarse-grained viscosity p* and heat conductivity k* are large enough 
in comparison to the microscale viscosity p and heat conductivity k (but at the same time 
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still small enough to enable the Navier-Stokes approximations in ( 5.9| >), one can simplify 
the above expression as 


(5.12) 


Ft=- 


Iv-^ 

2 P 


K 


*V ( ^ 


+ 2^* [V0M*] 


TS 


The main advantage of the coarse-grained Navier-Stokes equations in ( |5.10| > is that they 
are technically not very different from the usual Navier-Stokes equations, for which 
many computational methods have already been developed. Here, we have the addi¬ 
tional transport equation for the turbulent energy transport in ( |4.31[ ), but, again, it is 
not much different from the usual pressure transport equation, except that it has the 
additional forcing, dissipation, and diffusion. 

The main drawback of the coarse-grained Navier-Stokes approximation is, however, 
the following. Even if linear damping ( |5.7| > applies for a given type of turbulence (which 
does not restrict the use of the coarse-grained Grad equations in (5.61), the necessary 
condition for the coarse-grained Navier-Stokes approximation in (5.10) to be viable is 
that the evolution time scale of the coarse-grained stress and heat flux must be much 
faster than that for the density, velocity and pressure. This condition requires that both 
the coarse-grained viscosity and heat conductivity must be very small (which is usually 
the case for the conventional microscale }i and k). However, often in the models with the 
turbulent viscosity approximation these parameters can be several orders of magnitude 
larger than the microscale viscosity. In this situation, the time-scale separation is unlikely 
to exist, and thus the turbulent Navier-Stokes approximation may not apply. 

5.4. Coarse-grained Burnett and super-Burnett equations. A key advantage of the new 
coarse-grained variables in ( 3.5| ) is that they obey the well-studied hierarchy of the mo¬ 
ment equations in ( |3.6| ). Because of this, one has an opportunity to use existing clo¬ 
sure methods for these equations. Above we showed how to derive the coarse-grained 
Grad ( |5.6| ) and Navier-Stokes ( |5.10| ) closures for this hierarchy, but one does not have 
to stop there. Using the ratio }i* /p* (or, equivalently, k*/ p*) as a small parameter akin 
to the usual Knudsen number in the molecular kinetics, one can carry out the standard 
Chapman-Enskog perturbation expansion into the higher orders, obtaining the coarse¬ 
grained Burnett 0 and super-Burnett |32| equations in the same way they are derived 
for the conventional transport equations. While these equations are generally ill-posed, 
some relaxation-type regularizations were developed for them |T7|. Here we do not 
elaborate on this further, as the derivation itself is straightforward, while a more de¬ 
tailed study of the resulting coarse-grained Burnett or super-Burnett equations is an 
entirely separate topic. With the higher-order Chapman-Enskog expansion, the expres¬ 
sions for S* and q* in ( |5.9[ ) will contain additional terms, scaled by powers of p* /p* and 
K* /p*, which will accordingly affect the expression for in (5.11). The transport and 
dissipation parts of the equation for Ef in (4.31) will remain the same. 


6. Practical considerations for a computational implementation 

In this section we discuss possible issues arising in the course of numerical implemen¬ 
tation and modeling of the transport equations, developed above. 
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6.1. Inclusion of the coarse-grained spatial scale information. A key subject which was 
thus far left out of the picture is the spatial scale information about the coarse-grained 
and turbulent scales. While this scale information can be avoided in the formal coarse¬ 
grained transport formulation presented above (simply due to the fact that the Boltz- 
marm equation is linear except for its collision operator), it naturally resurfaces when 
one has to discretize the transport equations for a numerical simulation. For example, if 
the discretization scale is so fine that it resolves the viscous molecular dissipation, then, 
clearly, the turbulent energy Et must be (nearly) zero, as should be the coarse-grained 
stress S* and heat flux q*. On the other hand, if one chooses to coarsen the spatial 
discretization mesh, then Et, S* and q* should naturally become larger, which means 
that their corresponding dissipation terms must be weakened appropriately. The only 
way to set these conditions in the transport equations is to adjust the turbulent energy 
dissipation rate pet and damping in the coarse-grained collision terms Cg and c* (or, 
alternatively, coarse-grained viscosity p* and heat conductivity k*), solely because there 
is nothing else in the transport equations that can be changed. 

At first, this situation appears to conflict with the way the Reynolds average is initially 
set up as a statistical average, since the ensemble averaging is formally unrelated to the 
spatial filtering (like the one in the LES). To clarify the situation, here we offer an informal 
explanation how the spatial scaling information can be incorporated into the Reynolds 
averaging. The Reynolds averaging process consists of two stages - first, "ensemble 
generation", and, second, "ensemble averaging" (as purely a thought process, of course, 
since there is no explicit numerical ensemble generation or averaging involved). Here 
we claim that the spatial scale information is taken into account during the "ensemble 
generation" stage. Indeed, observe that all ensemble members must share the coarse¬ 
grained scale features (as it is the "resolved" scale), while, of course, differing on the 
turbulent scale. Therefore, if the coarse-grained scale is in fact so fine that it resolves even 
the molecular viscosity effects, the corresponding statistical ensemble must, essentially, 
consist of a single member (or many identical members), as any difference between the 
ensemble members on the unresolved scale is nearly instantaneously damped to zero 
by the molecular viscosity. Conversely, if the coarse-grained scale is not too fine, then it 
lifts the difference between the ensemble members into larger spatial scales, and hence 
weakens their dissipation. 

This leads to the natural conclusion that the spatial scale information must be encoded 
into the dissipative terms parameterizing the turbulent energy decay rate pet and the 
coarse-grained collision terms Cg and c*. Particularly, in the limit as the coarse-grained 
scale becomes fine enough to resolve the molecular viscosity scale, the turbulent energy 
decay rate pet must be strong enough to ensure that the turbulent energy Et is (almost) 
zero, while the coarse-grained collision terms Cg and must transform into their fine¬ 


grained analogs in (2.23). 


6.2. Numerical stability and possible shock formation. From was is developed above, 
one can envision three general scenarios of a practical computational set-up: 


(1) The coarse-grained Grad transport equations ( |5.6[ ) with the specific, problem- 
dependent parameterization of the collision terms; 
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(2) The coarse-grained Navier-Stokes transport equations ( 5.10| > with the coarse-grained 
viscosity and heat conductivity; 

(3) The coarse-grained Grad transport equations ( |5.6 t with collision terms param¬ 
eterized via linear damping as in (5.7), with coarse-grained viscosity and heat 
conductivity 

There is not much to comment on in the first two scenarios, because the latter is com¬ 
prised by the well-studied Navier-Stokes equations, while the computational proper¬ 
ties of the former heavily depend on the implementation of the coarse-grained collision 
terms. The third scenario is, however, of more interest, since it is typically used in rar¬ 
efied gas dynamics |T4|T5p^^^. In the turbulent flow simulations, the third scenario 


is likely to appear when the turbulence is still weak enough to parameterize the colli¬ 
sions by linear damping ( |5.7| ), but at the same time the coarse-grained viscosity and 
heat conductivity k* are not small enough to enable the Navier-Stokes approximations 
of S* and q* via ( |5.9[ ). In this situation, observe that the resulting transport equations for 
the coarse-grained variables do not have any diffusion terms, which leads to the creation 
of shocks at high Mach numbers |14|. 

In this scenario, one might consider the regularization strategy proposed by Struchtrup 
and Torrilhon |35|37|, which is, roughly put, a type of the Chapman-Enskog expansion 
built around the Grad state in ( 2.29| ), and applied to the transport equations for Q* and 
R*. The way this regularization works is the same as in the Navier-Stokes parameter- 
izations in ( |2.25| )- p.27| ), yielding additional diffusion terms in the transport equations 
for the coarse-grained stress ( |5.6d| ) and heat flux ( |5.6e| ), and thus dissipating the shock 
formation. The diffusion from these extra terms manifests on the scale fi* /p* (or, alter¬ 
natively, K* / p*), and extends to the larger scales as the coarse-grained viscosity and heat 
conductivity increase, thus allowing for coarser computational discretization. 


6.3. A crude approximation for the turbulent energy dissipation rate. The simplest 
way to model the turbulent energy dissipation rate pst is by using the dimensional anal¬ 
ysis. Above, we concluded that a parameterization for pst should depend on a length 
scale parameter which characterizes the coarse-grained spatial scale (such as the spatial 
discretization size). Let us now assume that et may only depend on such a spatial scale 
parameter L, and the turbulent energy Et itself, scaled by the density p: 

( 6 . 1 ) etoc(Et/pfL^ 


where the constants a and j6 must be determined by considering physical dimensions. 
Now, observe that the dimension of Et/p is length^/time^, while the dimension of et is 
length^/time^. Thus, the only way to equalize dimensions in the expression above is to 
set a; = 3/2, j6 = —1, which yields 


( 6 . 2 ) 


et oc 


r3/2 


^3/2L' 


pet 


cx 


r3/2 


p^/^E- 


This well-known simple parameterization was proposed by Taylor |36|. Observe that it 
is in agreement with what was concluded earlier; namely, it disappears when Et ^ 0 
(and thus does not let Et become negative), and increases for fixed Et as E ^ 0, thus 
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more strongly driving Et towards zero in the case the coarse-grained spatial resolution 
scale becomes more refined. 

However, we can immediately see that the simple parameterization above has a serious 
drawback. Recall that we agreed in Section |3 that the new coarse-grained transport 
equations must have solutions which do not exhibit rapid oscillations at small scales, 
otherwise the very point of introducing the coarse-grained equations would become 
moot. Nonetheless, with the turbulent energy decay rate pet given as above in ( |6.2| >, the 
turbulent energy transport equation in (4.311 has no diffusive terms of sufficient strength 
(the existing diffusion ferm there manifests on the molecular scale), and in that respect 
the equation in (4.311 is no different from the conventional transport equations in ( |2.7| > 
we started with. 

Therefore, here we contend that the turbulent energy decay rate pet must incorporate 
a diffusion term which manifests on the coarse-grained scale. In particular, we propose 
a simple extension of ( 6.2| > based, again, on a dimensional argument: 

, X > 

(6.3) pet oc 


—h ocL/^Et 


where a is a scalar non-dimensional parameter. Observe that above the first term is 
exactly the Taylor parameterization from \6.2) . The additional term, scaled by oc, is a 
diffusion term, which, contrary to the linear damping term, increases when L becomes 
large, extending the diffusion scale proportionally fo the coarse-grained scale parameter 
L. Conversely, the diffusion term disappears as the coarse-grained scale becomes more 
refined. This seems fo be in consistence with what one would expect from the behavior 
of the turbulent energy. 

Of course, the parameterization in ( |6.3| > is the simplest model for pet one could pos¬ 
sibly come up with, but at least it seems to be a reasonable starting point for use in 
computational modeling. There also exist more sophisticated models for et (such as, for 
example, the k-e and k-co models |38|) but they require a separate transport equation for 
et, which, in turn, introduces additional closure problems via unspecified parameters. 


6.4. Crude approximations for the coarse-grained collision terms. Here we offer a 
very crude parameterization of the coarse-grained collision terms in the form of a lin¬ 
ear damping, which, of course, will nof stand any critique when the properties of fhe 
coarse-grained collision terms become better studied, but at present may at least serve 
as a starting point in the computational modeling. Here we assume that the coarse¬ 
grained collision terms and c* can be modeled by a linear damping of the form ( |5.7[ >, 
where the coarse-grained viscosity p* and heat conductivity k* are to be approximated 
somehow. From the prior considerations, we expect that p* and k* must be proportional 
to the coarse-grained length scale (that is, they must become small when the length 
scale is short, and vice versa). At the same time, we know that, as the coarse-grained 
length scale becomes short, p* and k* must become their microscale counterparts p and 
K, respectively. The form of the microscale p and k is given in Grad p4| as 


(6.4) 


m 

P,K Oi — 



Pm^ 
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where m is the mass of the molecule, a is its linear size (diameter or radius), and pm 
is the "density" of a molecule. Now, with help of a vivid imagination, one could in¬ 
terpret the collision processes on a coarse-grained scale as ones with rather large and 
soft molecules. This means that for the coarse-grained viscosity or heat conductivity one 
likely has to keep the "density" pm fixed, but increase the "size" of the molecule to be 
that of the characteristic length of the coarse-grained scale, while at the same time also 
parameterizing p via either its coarse-grained counterpart, or its Reynolds average: 

' ,*\l/2 


(6.5) 


U*,K* OC ( ^ 

\P 


1/2 


L OC < 


-) 

p J 

pj 


or 


3p* - 2 Et 
3p 


1/2 


L, 


where the choice is between the lack of feedback from the small scales, or a feedback 
which damps the viscosity and heat conductivity if the turbulent energy becomes too 
large. The coefficients of proportionality should likely be chosen so that in the limit as 
the resolved scale becomes the molecular viscosity scale, both p* and k* become their 
microscale counterparts. In fact, they do not necessarily have to be constants, as that 
would impose the same Prandtl number on the coarse-grained scale, which may not 
necessarily be the case in practice. 

The relation above constitutes the simplest, crudest parameterization for the coarse¬ 
grained viscosity and heat conductivity. This parameterization is, however, consistent 
with what is expected of it, namely, it increases the coarse-grained viscosity and heat 
conductivity in proportion to the coarse-grained spatial scale, and, if the turbulent en¬ 
ergy feedback option is chosen, it also tends to decrease them as Et increases, thus pro¬ 
viding an additional balancing mechanism. Note, however, that a regularization along 
the lines of Struchtrup and Torrilhon |35|^ is needed for a linear damping parameter¬ 
ization in ( |6.5| > to improve the numerical stability of the coarse-grained Grad equations 
in ( |5.6| >. Alternatively, if the resulting p* and k* are found to be small enough to allow 
the steady-state parameterization of the kind in ( |5.9[ ), then one can optionally use the 
coarse-grained Navier-Stokes equations ( |5.10| > instead of the coarse-grained Grad equa¬ 
tions ( |5(6| > (or suitably regularized Burnett/super-Burnett expansions, for better accuracy 
with not-so-small p* and k*). 

7. Summary and future research 

In this work, we propose a new framework for the coarse-grained transport of a tur¬ 
bulent flow, which is based on the Reynolds averaging of the Boltzmaim equation in 
( 2.1|>, ra ther than the conventional Reynolds averaging of the Navier-Stokes equations 
in ( 2.28| >. The proposed framework consists of the equations in (3.6) for the transport 
of the coarse-grained variables of the flow in ( 3.5| ), and the equation for the transport 
of the turbulent energy in (4.31) under the approximation of the anelastic turbulence. 
We also propose two different closures for the hierarchy of the transport equations for 


the coarse-grained variables: the more general Grad closure in (5.6), and the more sim¬ 
ple Navier-Stokes closure formulation in ( 5.10| ), which can likely be expanded into the 
higher Burnett or super-Burnett orders to improve accuracy. The proposed transport 
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model includes three unknown parameters: the turbulent energy dissipation rate pst in 
the turbulent energy equation ( 4.31| >, the coarse-grained stress collision operator Cg (or, 
equivalently, the corresponding coarse-grained viscosity in the coarse-grained stress 
equation ( |5.6d| >, and the coarse-grained heat flux collision operator c* (or, equivalently, 
the corresponding coarse-grained heat conductivity k*) in the coarse-grained heat flux 
equation ( |5.6e| >. These three parameters caimot be defined exactly within the scope of 
the new transport model, as they depend on the properties of the nonlinear fluid-fluid 
interactions, and require separate treatment. We also suggest crude parameterizations 
for pet, Cg and c* based on dimensional analysis, to be used as a starting point in the 
practical computational modeling. 

Future research. The main deficiency of the developed coarse-grained transport frame¬ 
work is that it is derived under the assumption of a monatomic ideal gas model. While 
this assumption greatly simplifies calculations, at the same time it severely restricts the 
applicability of the framework, as most gases in the surrounding nature are at least di¬ 
atomic (e.g. the air). At the same time, such a framework would likely be of use in 
the applications which involve the circulation of the large scale atmosphere, as it is well 
known that the atmosphere is a highly turbulent medium, and at the same time present 
computational limitations restrict the spatial resolution of the global circulation models 
to rather coarse meshes. Thus, the development of the analogous coarse-grained trans¬ 
port framework for the polyatomic (or at least diatomic) ideal gases is presently our 
main priority. 
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